typeset using JPSJ.sty <ver.l.Ob> 



Dynamical Properties of the One-Dimensional Supersymmetric 
t-J Model: A View from Elementary Excitations 

Yasuhiro SaigaQ and Yoshio Kuramoto 1 

Institute for Solid State Physics, University of Tokyo, Roppongi 7-22-1, Minato-ku, Tokyo 106-8666 
1 Department of Physics, Tohoku University, Sendai 980-8578 

(Received ) 

Dynamical properties, such as dynamical spin and charge structure factors and single-particle 
spectral functions, are studied for the one-dimensional supersymmetric t-J model with inverse- 
square interaction. Exact diagonalization and the recursion method are used for finite systems up 
to 16 sites. A simple rule is proposed to understand the supermultiplet structure in the excitation 
spectrum. Numerical calculations show that the dynamical spin structure factor is independent 
of the electron density in the region where only two spinons contribute. In the electron removal 
spectrum from half-filling, close correspondence is found between the dominant contribution from 
three-spinon plus one-holon states in the t-J model and the two-spinon excitations in the Haldane- 
Shastry model. Comparison is made with dynamics of the nearest-neighbor supersymmetric t-J 
model. 

KEYWORDS: dynamical spin structure factor, dynamical charge structure factor, single-particle spectral function, 
t-J model, supersymmetry, inverse-square interaction, nearest-neighbor interaction, exact diagonal- 
ization, recursion method 



§1. Introduction 

One-dimensional (ID) electron systems show strong 
quantum fluctuations, and have been investigated inten- 
sively both theoretically and experimentally. Most of 
real quasi-lD systems have localized electrons and are 
regarded as spin systems. Recently some systems with 
doped holes have also been studied experimentally. For 
example, Yb4As3 has a small number of carriers, and 
shows a heavy-fermion behavior., in resistivity, magnetic 
susceptibility and specific heat. El The spin dynamics in 
Yb4As3 probed at low temperature suggests that the na- 
ture of low-energy spin excitations is similar, to that of 
the ID antiferromagnetic Heisenberg model. aa On the 
other hand, a compound Sri4Cu2404i consists of two-leg 
CU2O3 ladders and edge-sharing Cu02 chains. The num- 
ber of holes in the jphain is estimated to be 60% of the 
number of Cu ionsu-' One may ask how the existence of 
holes affects dynamical properties of quasi-lD systems. 
In the presence of both spin and charge degrees of free- 
dom, one can expect much richer physics as compared 
with spin-only systems. 

Recent angle-resolved photoemission measurements 
have revealed the existence of spin-charge separation, in 
quasi-lD systems. In Mott insulators such as SrCu02t3 , B' 
and Si'2Cu03,lP angle-resolved photoemission data show 
two distinct dispersions in the energy-momentum space. 
One is assigned to the holon band, and the other is re- 
garded as the spinon band. This is an indication of 
the spin-charge separation. More recently, an angle- 
resolved photoemission study on PrBa2Cu30y has been 
performed© Signals from the doped CUO3 chain re- 
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veal the spinon and holon bands in the range of mo- 
menta from to 7t/4 in units of the reciprocal lattice 
parameter. This indicates that the CUO3 chain is nearly 
quarter-filled. In order to investigate theoretically dy- 
namics of-these systems with holes, one can take the ID 
t-JmodelBllB' as the simplest model. To understand the 
photoemission data, in particular, stopWiXifJaakjirppa- 
gators in the t-J model is necessaxyMBByBBEP 

For the t-J model with nearest-neighbor interaction 
and small J/t, which is ordinarily adopted, dynamical 
correlation functions are rather intricate. This com- 
plexity comes from the fact that Landau quasi-particles 
are no longer well-defined elementary excitations. It is 
known that elementary excitations in ID t-J model obey 
fractional statistics. In this context, the long-range su- 
persymmetric t-J modelES' is suited for extracting the 
essential dynamical properties since the elementary ex- 
citations do not interact except for the generalized ex- 
clusion principleO-' For the long-range t-J model, the 
excitation contents are composed of spinons, holons and 
antiholons] the number of elementary excitations con- 
tributing to dynamics is finite due to the high symmetry 
of the model. Ef Ha and Haldane proposed regions of 
nonvanishing spectral weight in the energy-momentum 
space, which is called the "compact .support" , for vari- 
ous dynamical correlation functions.^ Recently, exact 
results far _thc spectral weight itself have been partially 
derived. EIr c3 However, full analytic solution of the dy- 
namical correlation function is not yet found. 

In this paper, we clarify the dynamical properties of 
the ID supersymmetric t-J model with long-range in- 
teraction. We calculate the dynamical spin and charge 
structure factors and the single-particle spectral func- 
tions for systems up to 16 sites with various hole num- 
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bers via exact diagonalization and the recursion method. 
Using the asymptotic Bethe-ansatz equations and the 
skew Young diagram, we reexamine—the excitation con- 
tents proposed by Ha and HaldaneHJ The relationship 
between the calculated dynamical quantities and the ex- 
citation contents is investigated. The results of dynam- 
ical quantities are compared with thx^e for the nearest- 
neighbor supersymmetric t-J modelj 2 ^ which is solvable 
by the Bethe ansatz. 

The remainder of this paper is arranged as follows. 
In §2 we consider two boundary conditions for the 
long-range t-J model: the periodic boundary condition 
and the antiperiodic one. Particularly, the antiperiodic 
boundary condition is nontrivial for the long-range t-J 
model. In §3 we show results of various dynamical quan- 
tities and interpret them in terms of elementary excita- 
tions. Finally we give the summary of this study and 
discuss the results in §4. 

§2. Model and Boundary Conditions 

We consider the ID t-J model 



JV 



"=£ E 



2<j £ — — oo 



-tijt (4*Cj+eN,* + H.c. 



+Jije [ Si ■ Sj+tN — -riiiij+eN 



(2.1) 



with even number N of sites and average electron number 
h per site. Here c ia = c ia {\ — n it _ a ) is the annihilation 
operator of an electron with spin a at site i with the 
constraint of no double occupation. The summation with 
respect to I has ibeea introduced to treat the system in 
a ring geometry.EjEs' Two types of the interaction are 
considered: 

(i) the nearest-neighbor (NN) type, i.e., tijt = tSj^+iSe^; 

(ii) the inverse-square (IS) type, i.e., tijt = t/(xi — Xj — 
IN) 2 . 

The Hamiltonian in the cases (i) and (ii) are written as 
7Ynn an d His, respectively. 

For the NN t- J model, it is conventional to choose the 
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periodic boundary condition Cj 

2 with iV c being the electron number and m an integer 
and the antiperiodic boundary condition Cj+N,a = —Cjo 
for N e — 4m. This leads to the nondegenerate ground 
state which is a singlet with zero total momentum. The 
Hamiltonian is written out as follows: 



JV 



-*E (C^+i 



H.c. 



-J I Si' Si+i - -mrii+i 



(2.2) 



where cjv+i )CT = c\ a for iV c = 4m + 2, and cn+i.o- = —c\o 
for N e = 4m; in addition, Sn+i = Si and n^ + i = ri\ in 
any case. 

For the IS t-J model, on the other hand, the effect 
of the boundary condition is nontrivial. When the pe- 



riodic boundary condition Cj+tN,o 
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are imposed, 



both transfer energy and exchange energy become the 
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tf=tD{ Xi 



(2.3) 
(2.4) 



Under the 



where D(xi — Xj) = {N/tt) sin[7r(a;i — Xj)/N 
antiperiodic boundary condition Cj-^tN a- = {—L)"Cj a , on 

(a) 

the other hand, the transfer energy t\j is given by 



,(a) 



=*E 



£N) 2 



= tD(xi - x i )" 2 cos[7r(x i - Xj)/N]. (2.5) 

This corresponds to., the special case of the twisted 
boundary conditionEs- 13 

We want to consider the spin and charge dynam- 
ics from the nondegenerate ground state. Hence for 
N e = 4m + 2 we choose the periodic boundary condi- 
tion, and for N c — 4m we impose the antiperiodic one. 
Consequently, we consider the following Hamiltonian as 
the t-J model with IS interaction: 
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-Ji: 
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where is given by eq. (2.4), and 
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ij 

,(a) 



for N c = 4m 
for N c = 4m. 
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When we calculate the single-particle spectral functions 
from a less-than-half-filled ground state, we start from 
the initial state with N c = 4m + 2 and therefore use 
Uj = t:J' . We always take the supersymmetric case t = 
J/2 in the following calculation. 

The dynamical quantities we consider in this paper 
are the followings: the dynamical spin structure factor 
S(q,u), the dynamical charge structure factor N(q,u), 
and the single-particle spectral function A(k,u>). These 
quantities are defined as follows: 



S"(q,u) 
N{q,uj) 
A-(h,w) 



'^\(u\a'JO)\ 2 8(u-E v + Eo) ) (2. 



' S(lj - E v + E ), (2.9) 
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N-^^sfe-^, n q = N-Vi^im- 



i)e"^, and c ka = A^ 1/2 £, 



ikt 



; q (or k) and u 



represent the momentum transfer and the excitation en- 
ergy, respectively. \v) denotes an eigenstate of H with 
energy E v with Eq being the ground-state energy. These 
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dynamical quantities can be written in the form of a 
continued fractionBS' We truncate the continued frac- 
tion after 100-200 iterations. In our calculation, poles 
and intensities have accuracies of 7 digits and 4-5 digits, 
respectivelyr$ 

§3. Dynamical Properties and Elementary Ex- 
citations 

3.1 Ha-Haldane's support 

Before presenting our results on dynamical proper- 
ties, we briefly review the compact supports for vari- 
ous correlation functions for the IS supersymmetric t-J 
model. The compact support means a finite region where 
the spectral function is nonzero. To draw the support, 
one needs three pieces of information:E2P dispersion re- 
lations of related elementary excitations, excitation con- 
tents of the intermediate states for each correlation func- 
tion, and selection rules. The dispersion relations for the 
right (left)-going spinon (sr/l)), holon (h R (L)) and anti- 
holon (h) in the thermodynamic limit are given by 



£ SR(L) A : 

e hR(L) A 

1 



(3.1) 
(3.2) 

(3.3) 



where v® — ir and v® = n(l — n). The range of mo- 
mentum permissible for each elementary excitation de- 
pends on the electron density n. The right (left) spinons 
and holons are allowed in < p < kp (— fcp < p < 0), 
where fcp = nn/2. The antiholons propagate in the re- 
gion — v® < p < v®. The excitation contents, which 
mean all possible states excited from the ground state 
by a local operator, are determined with the help of the 
numerical results. Finally, the selection rules are found 
empirically, also with use of the numerical information. 

3. 2 Derivation of excitation contents 

Although exact expression of the matrix element, or 
the form factor, is not yet obtained, one can determine 
the excitation contents by comparing numerically ob- 
tained levels with those in the asymptotic Bethe-ansatz 
equations. In this subsection, we explain how to de- 
rive the excitation contents corresponding to each energy 
level. 

Let us take the completely up-polarized state as the 
reference state and consider M + Q pseudo-particles 
which represent M electrons with spin down and Q holes. 
Note that these pseudo-particles are not spinons and 
holons. The aspaap-totic Bethe ansatz leads to the fol- 
lowing equation:E3B 



E 
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-Q l 
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where e(p) = p(p — 2n)/2. The pseudomomenta p^ are 
determined by the following equations: 

Q 

p^N = 2tt J,, - it J~] sgn(p M - qi) 



M+Q 

+7T ^2 s S n (Pn-Pv), 

v=\ 

M+Q 

27rl; = 7T ^ s S n (* -Pli), 
H=l 



(3.5) 
(3.6) 



where J M £ \{M + l)/2,N - (M + l)/2] for \i = 
1, 2, • • • ,M + Q, and h £ [-{M + Q)/2, (M + Q)/2] 
for i — 1,2, ■ ■ ■ ,Q. The quantum numbers and Ii 
are integers or half-integers, and are arranged in the as- 
cending order. We assume that all p^ and qi are dis- 
tinct from each other. For given p^, we specify such 
that < Pp. < qj(ji)+i in the case of < < Q, 
j(p) = if p^ < qi and j(p) = Q if q Q < p^. Then we 
obtain 



N 
2^ 



Pa* = j m 



AI + 1 



-iO)> 



(3.7) 



from eqs. (3.5) and (3.6). Note that p^ = Np^/^n) is 
always integral. 

Identification of excitation contents proceeds as fol- 
lows: First, we give arbitrary configurations of J ^ and 
Ii in the subspace with fixed M and Q and obtain 
a set of pseudomomenta {p^}. Second, using the set 
{p^}, we find the total energy E, the total momentum 
P = Y^!=\^ P^ ( m °d 2tt), and the sequence of occu- 
pied state (represented by 1) and empty state (repre- 
sented by 0) in the pseudomomentum space. This se- 
quence is referred to as the motif. The first and last 
entries of the motif are always set to 0. We note that 
e{Pfj) takes the maximum value (= 0) for p^ = and 
2n. Finally, the sequence (i.e., the motif) is assigned 
to an eigenstate with the same energy and momentum 
obtained by exact diagonalization. Let us give an exam- 
ple. For (N,Q,M) — (4,2,0), permissible values of J M 
are 1/2, 3/2, 5/2 and 7/2, while those of h are -1,0 
and 1. We consider the case of {J M } = {3/2,7/2} and 
{Ii} = {—1,0}. Then we obtain {p^} = {1,3}, so that 
the motif becomes 01010. 

The IS supersymmetric t-J model has the supermul- 
tiplet structure, ija.. enormous degeneracies in the ex- 
citation spectrumEiP However, the motif itself does not 
provide us information on the degeneracy. To extract the 
information, the skew Young diagram methodEj is very 
useful. The motif and the skew Young diagram have a 
one-to-one correspondence. The correspondence goes as 
follows: To begin with, we write the first box. Then we 
read a motif from the left, neglecting the first and last 
entries of the motif; if we encounter 1 (0), we add a box 
over (on the right of) the last box. We now propose a 
simple rule to reproduce the supermultiplet structure of 
the IS t-J model by, generalizing the known rule for the 
SU(m) symmetry.E9 Namely one puts an index for the 
internal degrees of freedom on each box of a skew Young 
diagram. In the present case we represent an up-spin by 
1, a down-spin by 2, and a hole by o. The order of them 
is denned as 1 < o < 2. We introduce the following rule: 

[i] a < b if b is lower-adjacent to a; however, consecutive 
o's are allowed in a column. Note that the numbers 1 
and 2 can appear at most once in a column. 

[ii] a > b if b is left-adjacent to a; however, consecutive 
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o's are forbidden in a row. 

For the 4-site IS t-J model with the fixed number of 
holes, the degeneracy based on this rule is consistent 
with the numerical result for all eigenvalues. We note 
that this rule enables us to decompose the degeneracy 
for each number of holes. For instance, the motif 01010 
for N = 4 has the energy E/t = ir 2 (5Q - 12)/16 and 
the momentum P = 0. After translating this motif into 
the skew Young diagram, we find that the number of 
ways to put 1, 2 or o on each box is 12, according to 
[i] and [ii] above (see Fig. 1). The 12-fold degeneracy is 
decomposed as follows: nondegenerate for Q = 0, 4-fold 
degenerate for Q = 1, 5-fold for Q — 2, and 2- fold for 
Q = 3. These results agree with those via numerical 
diagonalization: (Stot = 0) 1 for Q — 0, (Stot = 1/2) 2 
for Q = 1, (S tot = 0) 2 9 (Stat = I) 1 for Q = 2, and 
(Stot = 1/2) 1 for Q = 3. 

We note that the degeneracy in the spectrum has also 
been studied for the SU(m)EJ and SU(m|n) supersym- 
metrico* Polychronakos spin chains, which are a variant 
of the l/r 2 -type model. The SU(1|2) specialization of 
the latter chain has the same supersymmetric Yangian 
symmetry as the IS t-J model. Thus one may expect 
that the supermultiplet structures in the spectrum are 
the same for both models. We have checked the follow- 
ing for N — 4: If the first term of eq. (3.4) is subtracted 
from an energy level by exact diagonalization of the t- 
J model, the degeneracy of the energy level is consistent 
with that for the SU(1|2) supersymmetric Polychronakos 
spin chain.B The subtraction is necessary since the first 



term of eq. (3.4) breaks the global SU(1|2) supersymme- 
try. 

Now, we can extract the excitation contents from 
the motif. According to Ha and HaldaneJ 2 !! a spinon 
(with spin 1/2 and no charge) is represented by two 
consecutive zeros, and a holon (with charge — e and no 
spin) by two consecutive ones. For less-than-half-filling, 
the motif of the ground state with no magnetization 
is 01010. ..10|111...111|01... 01010, which has (N + 1) se- 
quence. The central region [(Q + 1) consecutive ones] is 
regarded as a "pseudo Fermi sea (PFS)" . A hole in the 
PFS (i.e., zero in the center region) is called an antiholon 
(with charge +2e and no spin). Spinons and holons are 
excited in the left and/or right region. For example, a 
motif 00110|1111101|10100 should read (s L , h L ) + h + 

(hR, SR.). 

When we derive the excitation contents in this 
way, however, exceptional situations occur in some 
cases. For motifs such as 01010|1111101|01100 and 
01010|1111110|11010, which are actually present for 16 
sites with 6 holes, excitation contents according to ref. 
20 become h + (Iir, sr) [case (a)] and h + Iir [case 
(b)], respectively. These contents do not fulfill the 
charge and/or spin conservations, and the latter does 
not fulfill the semionic statistics of holons. This sug- 
gests that the excitation contents by Ha and HaldaneEJ 
do not always hold true for finite systems. In contrast, 
without.. holes, i.e., for the Haldane-Shastry (HS) spin 
model J2j'l3 we can always read two spinons from the 
motif for all states that contribute to S(q, ui) even for 
finite systems (N < 16). 



Excitation contents for a finite system can be extracted 
also by the skew Young diagram. We first determine the 
PFS of the diagram from the corresponding motif. Note 
that the PFS consists of Q boxes. We regard o outside 
the PFS as a holon, and a row of a 1-2 pair in the PFS 
as an antiholon. Single 1 or 2 in a column, both of which 
do not make a pair, is regarded as a spinon; it does not 
matter whether the column includes o or not. Then, in 
the case (a), some of the diagrams with the entries of 
1, 2 and o can read (sl, Iil) + h + (hR, sr). One of 
them is shown in Fig. 2(a). In the case (b), on the other 
hand, we find from Fig. 2(b) that only 1 of a 1-2 pair in 
a row is present in the PFS. We shall call such a part 
of the 1-2 pair in the PFS a "half-antiholon (h*)". This 
quasi-particle is assumed to have charge +e, no spin and 
semionic statistics. Thereby the diagram in the case (b) 
reads h* + hR, which satisfies both charge and spin con- 
servations and the statistical feature. In the case where 
an antiholon in the Ha-Haldane's excitation contents is 
replaced by a half-antiholon, a holon should be simulta- 
neously subtracted from the contents. As the number of 
sites increases with hole density kept constant, the effect 
of the edges in the PFS should become negligible. This 
means that a half-antiholon has zero energy. In fact, 
for N(q,u>), the energy levels with contents h* + hR are 
almost along the holon dispersion for q < k-p in the com- 
pact support. Thus only Ha-Haldane's contents should 
survive in the thermodynamic limit. 

We notice that the spinons and holons identified above 
do not satisfy the statistical matrix deduced from Ji- mi- 



croscopic derivation of the exact thermodynamicsEiP at 
least for N — 4. Nevertheless, these quasi-particles obey 
the semionic statistics and exhaust the Hilbert space for 
the t-J model. This reflects the fact that a different set 
of quasi-particle picture is available in spanning the com- 
plete Hilbert space. 

For a dynamical correlation function of the finite-size 
IS t- J model, an energy level with the Ha-Haldane's exci- 
tation contents has not always finite intensity. For exam- 
ple, let us take an eigenvalue with (q,u>) — (it, 1.09662t) 
for (N, Q, M) = (12,2,5). This eigenvalue has 4-fold 
degeneracy: (Stot = 0) 1 © (Stot = I) 1 , an d its motif 
00101|101|10100 leads to (s L , h L ) + h + (h R , s R ). Nev- 
ertheless, there is finite intensity for S(q,uj) but no in- 
tensity for N(q,ui). This implies existence of additional 
selection rules for N(q,ui). 

3.3 S(q,uj): dynamical spin structure factor 

We address to the spectral weight itself for various 
correlation functions in the remaining part of this sec- 
tion. Let us begin with the results of S(q,u>) for the 
IS t-J model. Figures 3(a) and 3(b) show the results 
in the 16-site chain for electron densities n — 0.875 
(14 electrons) and 0.5 (8 electrons), respectively. The 
solid lines show dispersion relations of elementary exci- 
tations for the compact support. These figures reveal 
clearly the compactness of the support for the excita- 
tion at arbitrary filling; there is no intensity outside the 
outermost dispersion lines. For two holes, the main in- 
tensity is in the low-energy region. As the hole den- 
sity becomes larger, the peak frequency at each q seems 



Dynamical Properties of the One-Dimensional Supersymmetric t-J Model 



5 



to approach the antiholon dispersion which runs over 
kp < q < ir. For n — 0.5 in Fig. 3(b), poles with domi- 
nant intensity are nearly along the antiholon dispersion 
over 0.25 < q/n < 1. 

Let us look at poles and intensities more carefully. A 
surprising feature is that poles and intensities in the two- 
spinon region agree with those for the HS model within 
the numerical accuracy. We have confirmed this fact for 
N < 16 with various n. It leads to the following conjec- 
ture: independent of n, the exact expression of S(q, lo) 
in the 2sB t -qnd 2sl regions is identical with that for the 
HS model^} i.e., 

S(q,Lo) =l/{A[(u- UL-(q))(u> - c L+ (q))] 1/2 } , 

(3.8) 

where wl-(<z) = tq{i: — q) and WL+(g) = t(q — 7r)(27r — q). 
This strong spin-charge separation in S(q, lo) should be a 
manifestation of the supersymmetric Yangian symmetry 
of the present model. Note that in the NN case such 
strong separation does not occur. 

In the region including holons and antiholons, on the 
other hand, the w-dependence of S(q, lo) is more compli- 
cated. Possible reasons for this complication are: 

• For a finite system, the intermediate states of S(q, lo) 
are not always (sl, ^l) + h + (ha, sr) states, in 
contrast with t he claim by Ha and HaldaneEJ As 
mentioned in §3.2, the intermediate states include 
the states with contents (sl, Iil) + h* + sr and the 
mirror ones (L R). 

• It is likely that the intensities of states with different 
excitation contents merge to a single pole. 

Now we consider the static spin structure factor 5(g), 
which is obtained by integration of 5(g, lo) over lo. Inde- 
pendent of the IS t-J model, Gebhard and VollhardtS 
calculated, S(q) associated with the Gutzwiller wave 
functionta' as follows: 



S(q) 



t 



i-I 



ln(l - n), 



for < q < 2k F , 
for 2fc F < q < 1"- 



(3.9) 

It was then demonstrated that the ground-state wave 
function of the IS t-J model is identical to the Gutzwiller 
oneJHf Thus eq. (3.E) turns out to be the exact expres- 
sion of S(q) for the IS t-J model. For < q < kp, 
integration of cq. (3.8) over lo reproduces the former of 
cq. ( |3.9| ). This supports our conjecture on S(q,Lo) in the 
two-spinon region . In Fig. 4 we show the results of 5(g) 
for N = 16 for various fillings. 

To compare with results for the IS model, we show 
the counterparts for the NN model. Figures 5(a) and 
5(b) show the results of S(q,u>) for two holes and eight 
holes in the 16-site chain, respectively. For two holes, 
the global features are similar to those for the IS model. 
However, there are two remarkable differences: 

[A] Contribution of higher-order spinons, holons and 
antiholons to S(q, lo) spreads over the high-energy region. 

[B] Even below the spinon dispersion, there are appre- 
ciable intensities for < q < kp. This corresponds to 
the lower holon dispersion in the NN model.Ef 

For eight holes, i.e., quarter-filling [see Fig. 5(b)], the 



region with main intensity in the (g, cj)-plane is rather 
different from the region with finite intensity for the IS 
type. This is attributed to the change of elementary exci- 
tations with increasing hole density in the NN model. c3 

3.4 N(q,Lo): dynamical charge structure factor 

We discuss the results of N(q,Lo) for the IS t-J model. 
Figures 6(a) and 6(b) show the results for 16 sites and 
various electron densities. All poles with finite intensity 
are inside the corresponding support which is indicated 
with the solid lines in the figures. For two holes in Fig. 
6(a), a typical example of lightly hole-doped case, the 
following features are observed: 

• The peak frequency of N(q, lo) for each q is at the 
upper edge of the continuum. 

• The intensity of lowest-energy at q — 2ir — 4fc F = 
7r/4 is stronger than that at q = 2kp = 7ir/8 for 
n = 0.875. 

As the hole density is increased, these features disappear 
gradually. 

According to Ha and Haldane,c3 the intermediate 
states of N(q, lo) consist of states with contents h + 21ir, 
the mirror ones, and the ones with contents (sl, he) + h 
+ (Iir, sr). We look at the intensity in the pure h + 21ir 
region, which does not overlap with the region including 
spinons. It is found that the intensity in the region in- 
creases monotonically with lo. We note that N(q,u>) in 
the region q < kp has recently been computed analyti- 
cally. EP 

Next, we discuss the static charge structure factor 
N(q). Like 5(g), the exact expression of N(q) is derived 
by Gebhacd and VollhardtxS' We define the following 
functions:c3 

at t \ x (1-n + x) 

Ni(x) = x - - In- 



(1 



N 2 {x) 



x) 



2 
ln(l 



(1-n) 



n 



x 



N 3 (x) = 2 - 2n 
N A = 2n + ln(l 
N 5 (x) = 2 

+ ln(l 



x (n — 1 + x) 
— In ■ 
2 



(1 



n 



x y 



n), 



2n 



l + x) 



x (n 

— In r- 

2 (1 + n-x) 



n 



x). 



(3.10) 

(3.11) 
(3.12) 
(3.13) 

(3.14) 



With use of these equations, when < n < 1/2, one has 

Ni(q/ir), forO<g<2fc F , 
N(q) = { A 2 (g/7r), for 2k F < q < 4fc F , (3.15) 
N 4 , for 4fc F < q <tt. 



For 1/2 < n < 2/3, one has 
N(q) = 



N^q/ir), for < q < 2k F , 
N 2 (q/n), for 2fc F < q < 2tt - 4fc F , 
N 5 (q/ii), for 2ir - 4fc F < q < it, 



(3.16) 
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and when 2/3 < n < 1, 



N(q) = { N 3 (q/n), 
N 5 (q/n), 



for < q < 27r-4fc F , 
for 2n - 4fc F < q < 2k F , 
for 2kp < q < tt. 

(3.17) 

It is clear that for n = 1/2 and 2/3 there are two different 
momentum regimes, otherwise there are three regimes. 
For < n < 1/2 and 4fcp < q < tt, we find the following 
relation between N(q) and S(q): 



N(q) + 4S(g) = 2n. 



(3.18) 



The sum N(q) + 4S(q) represents the_density correla- 
tion between the same species of spinier The fact that 
this quantity coincides with the noninteracting value 2n 
seems to be related to a hidden symmetry of the model. 
Figure 7 shows the numerical results of N(q) for 16 sites 
and various fillings. 

Finally, we compare the results of N(q, uj) for the NN 
model with those for the IS model. Figures 8(a) and 8(b) 
show the results for the 16-site NN t-J model with two 
holes and eight holes, respectively. For two holes, charge 
dynamics for the NN interaction shares the global feature 
of the IS case. However, there is a remarkable difference 
in contribution of higher-order spinons, holons, and anti- 
holons. This corresponds to [A] mentioned for S(q,uj) in 
the previous subsection. For eight holes, the excitation 
spectrum with main intensity deviates strongly from the 
compact support for the IS interaction—. This is due to 
modification of elementary excitations.El 

3.5 A(k,uj): electron removal from half-filling 

We examine the electron removal spectrum A~ (k, uj) 
from half-filling. For the IS t-J model, excitation con- 
tents consist of one-spinon plus one-holon (lslh) states 
and three-spinon plusj— one-holon (3slh) states in the 
absence of antiholons.EI Figure 9 shows the result of 
A~(k,u) for 16 sites. Note that uj = E (N C )-E U (N C -1) 
is negative [see eq. (2.10)]. The origin of energy is 
taken to be the ground-state energy at half-filling, and 
therefore the lowest energy at k — tt/2 is equal to the 
chemical potential — Eq(N c ) — Eq(N c — 1) for the 
finite-size system. We notice that /i~ and fj, + given by 
Eo(N c + 1) — Eq(N c ) are discontinuous in the gapped 
system; in the present case (i.e., Mott insulators with 
infinite energy for double occupation) /i + cannot be de- 
fined. The region with finite intensity is compared with 
the support for the one-particle Green function in the 
high-density limit n — > 1 (the solid lines of Fig. 9). In the 
thermodynamic limit, fj, (= dEo/dN e ) becomes — 7r 2 i/12 
(the broken line of Fig. 9). This value is consistent with 
— 0.822467i extrapolated from the numerical results for 
N = 4, 8, 12 and 16. The excitation at k = has the 
energy — 7r 2 i/3, which agrees with the extrapolated value 
(— 3.28987i) from the numerical calculation for N = 4- 
16. 

Let us turn now to the magnitude of the intensity. 
Recently, Kato has obtained tha.exact expression of the 
lslh contribution for A~(k, uj)E3c5> For a finite system, 
the lslh contribution Gi s ih(x,t) to the hole propagator 



at half-filling is given b"' 4 
1 

27V 

- N/2-2 N/2-K-2 



A(-iy 



T[N/2 - K+ 1/2}T[K + 1/2] 



ttN 2 ^ 

K=0 

-HE 



1=0 



T[N/2- K]T[K+ 1] 



xe -^iK cos (2 7ra; (/ + K + 1)/JV), 



(3.19) 



where E IK = 2tt 2 (2/ 2 - 21 M - 2K 2 + 2KM + I + K- 
M - l/6)/iV 2 + tt 2 /3 with M = N/2 - 1 (K and / 
being integers), and F[- • •] is the gamma function. From 
this result, we can read the momentum transfer fc, the 
excitation energy u> from the half-filled ground state, and 
the form factor Mi s ih as follows: 

k = 2tt{I + K + l)/iV-7r, (mod27r), (3.20) 

w = E IK , (3.21) 

2 T[N/2 - K+ 1/2}T[K+ 1/2] 



lh 



-. (3.22) 



irN T[N/2- K}T[K + 1] 

The lslh contribution for N — 16 is shown in Fig. 10(a). 
It is found thai at k = tt, the lslh state does not con- 
tribute at allS 

Analytical expression of the 3s lh contribution is not 
yet derived. We can derive numerically the 3slh contri- 
bution by subtracting the lslh contribution [Fig. 10(a)] 
from the total spectral weight [Fig. 9]. We show the 
3slh contribution for TV = 16 in Fig. 10(b). Note that 
there are some cases where the lslh and 3slh contri- 
butions merge into a single pole; for example, (k, lo) = 
(tt/4, -3.43123t) and (3tt/8, -4.89625i). 

Figure 10(b) gives us the following information about 
the 3slh contribution: 

• The 3slh contribution is present even in the region 
of the (k, tt>)-plane where the lslh states dominate. 

• The 3slh contribution has the peak intensity for 
each k near the lower boundary of the lslh con- 
tribution. This lower boundary for < k < n cor- 
responds to the spinon dispersion at half-filling. 

For each momentum, the pole positions of 3slh states 
with dominant intensity [the centers of black circles in 
Fig. 10(b)] find correspondence with the two-spinon exci- 
tations in the HS model. Namely, the values of — (uj — ujq) 
agree with the pole positions of S(q, ui) for the HS model, 
where uiq denotes the position of the single pole at k = 0. 
However, those intensities are always weaker than the 
counterparts for the HS model. These facts mean the fol- 
lowing: The 3slh state with dominant intensity consists 
of two moving spinons together with another spinon with 
zero energy and a holon with the maximal energy. The 
intensity is suppressed as compared with the two-spinon 
intensity, namely S(q, uj) in the HS model, because of 
addition of a further spinon and a holon. Nevertheless 
the rough feature of the intensity distribution resembles 
S(q,uj) in the HS model. 

Next we study the spectrum of electron removal in the 
NN t- J model, for which the result is shown in Fig. 11. 
Like the IS model, the number of poles for k — is only 
one. The remarkable differences are as follows: 

[a] There are many poles with small intensity outside 
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the 3slh continuum, although part of those poles are not 
apparent in the figure. 

[b] The branch of the lowest energy for < k < ir is 
the cosine-like function, in contrast with the quadratic 
function for the IS model. 

For the electron removal from half-filling, the elemen- 
tary excitation picture on the basis of spinons and holons 
applies only for the low-energy region in the NN super- 
symmetric t-J modelMd^ The difference between the NN 
case and the IS case appears in the dispersion relations 
of spinons and holons. Taking account of the spin and 
charge conservations, possible excitations for A~(k,to) 
consist of the _/V sp -spinon plus one-holon states where 
A sp = 1, 3, 5, • • • , N - 1. In the IS case the N sp > 3 
contribution vanishes. 

3.6 A(k,u>): single-particle spectral function at less- 
than-half- filling 
As the single-particle excitation from the ground state 
for n < 1, we consider the electron removal spectrum 
A~(k,uj) and the electron addition one A + (k,uj). The 
results for 16 sites with 14 electrons and 6 electrons 
are shown in Figs. 12(a) and 12(b), respectively. The 
chemical potential for the finite-size system is defined as 
yT = Eo(N c ) — Eq(N c — 1) for the electron removal, and 
as /i + = Eo(N c + 1) — Eq(N b ) for the electron addition. 
In the thermodynamic limit, /i - and /i + should be equal 
to jj, given by 



fl/t 



— 7r (3n 
12 V 



6n + 4). 



(3.23) 



This is .obtained by differentiating the ground-state en- 
ergy Ej^S> with respect to N e . The broken line of Fig. 12 
indicates the value given by eq. (3.23). The support by 
Ha and Haldane is shown by the solid lines. Our result 
for N — 16 reproduces their support. 

As the hole density becomes larger, the main inten- 
sity accumulates along the edge of the compact support 
and across the chemical potential line. This curve cor- 
responds approximately to the energy band in the low- 
density limit n — > 0. The intensity of A + (k,u) increases 
with increasing w at each k. On the other hand, the 
w-dependence of A~(k,Lu) is more intricate. The com- 
plexity arises because more elementary excitations con- 
tribute to A~(k,uj). Namely, the excitation contents for 
A + (k,oj) consist of (sl, Iil) + h and the mirror states, 
while (sl, Iil) + h + 2Xsr, hrt) and the mirror states 
contribute to A~(k,uj)!E$ 

Now we consider the sum rule. Because double oc- 
cupation at each site is excluded in the t-J model, the 
following sum rule holds: 



du [A + (k,w)+A-(k,u))] 



(3.24) 



independent of k (see Appendix). Integration of 
A^(k 1 uj) over lu yields the momentum distribution for 
spin cr, which can be defined as, 

n a {k) = (0\£{ a c kt7 \0), (3.25) 

where |0) denotes the ground state. For the IS t- J model, 
the exact expression of n a (k) is given by an infinite se- 



ries 

Before closing this subsection, let us compare with the 
single-particle spectral functions of the NN t-J model. In 
Figs. 13(a) and 13(b) we show the results of A ± (k, u>) for 
16 sites with two holes and ten holes, respectively. The 
prominent difference from the IS case is the curvature of 
the branch with dominant intensity representing cither 
a particle or a hole in the band picture. Again, we find 
contributions of higher-order spinons, holons and anti- 
holons in the high-energy region in the numerical data. 
For A ± (k,uj) with two holes, the area of appreciable in- 
tensity is similar to the compact support in the corre- 
sponding IS case. In A~{k,ui) with ten holes, the region 
withernam intensity shifts closer to the chemical poten- 
tialciP than that for the IS model. 

§4. Summary and Discussion 

We have investigated the dynamical properties of the 
ID supersymmetric t-J model for the two types of the 
interaction. The main results are summarized as follows: 

(1) For the inverse-square type, the dynamical spin 
structure factor is independent of the electron density in 
the spectral region where only two spinons contribute. 

(2) Concerning the electron removal spectrum from 
half-filling for the inverse-square type, the three-spinon 
plus one-holon states with dominant intensity have 
the correspondence with two-spinon excitations in the 
Haldane-Shastry model. 

(3) For the nearest-neighbor type with lightly doped 
holes, the region with dominant intensity for a dynamical 
correlation function is similar to the corresponding com- 
pact support of the inverse-square type. However, the 
regions for the two types become different as the hole 
density is increased. 

The feature (1) is an indication of the strong spin- 
charge separation in the zero-magnetic-field dynamics for 
the long-range supersymmetric t-J model. In applied 



field, the strong spir 
the charge dynamicsl 22 



irge separation appears also in 
i cB We will report detailed results 
on the dynamics in nonzero magnetic field in a separate 
paper. 

The feature (3) suggests the following conjecture about 
the dynamics for the nearest-neighbor type: In the ther- 
modynamic limit, the most relevant states for a correla- 
tion function would be composed of the excitatipa con- 
tents (i.e., the minimal sets) by Ha and HaldanejEI pro- 
vided the electron density is near half-filling. For large 
hole density, however, this does not hold. For the inverse- 
square type, the dynamics is characterized by spinons, 
holons and antiholons at arbitrary filling. On the other 
hand, for the nearest-neighbor type, proper picture of 
elementary excitations changes gradually with increas- 
ing hole density. Namely the Ha-Haldane's excitation 
contents no longer apply even approximately to the case 
with a large hole density. This is because the spinons 
and holons obeying the fractional statistics interact with 
each other. The charge carried by the holon and anti- 
holon excitations in the nearest-neighbor case is rather 
different from_that in the inverse-square case, even for 
quarter filling. O 
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Appendix: Sum Rule of A(k,uj) 

In this appendix we prove eq. ( |3.24[ ). Integration of 
A + {k,u) + A~(k, oj) over lu reduces to the quantity 



\ckac\o 



■4<7Cfc<7|0). 



(AT) 



Note that the tilde implies exclusion of double occupa- 
tion at each site. The Fourier transform: 
1 x - 



Cko 



leads to 



N ^ 

3 



cja = Cj>(1 - rij-v), (A-2) 



'j'a C 3<j) 

x(l-n i ,_ <r )(l -ny 



ikU'-i) 



= 1 - 



(A3) 



Here we have used the anticommutation relation of a 
fermionic operator Cj a . The expectation value of the 
second term of eq. (A-3) by |0) becomes n/2 since the 
numbers of electrons with a and —a are equal in the 
singlet ground state. 

This proof relies only on the condition that double 



occupation is excluded at each site. Thus eq. (3.24) holds 
irrespective of the interaction range and the value of J/t 
in the t-J model. 
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Fig. 1. Enumeration of all possible ways to put an index on each 
box of the skew Young diagram corresponding to the motif 01010. 
The entries of the numbers 1, 2 and o obey the rules [i] and [ii]. 



Fig. 9. A(k,w) of the IS supersymmetric t- J model with 16 sites 
at half-filling. The intensity of each pole is proportional to the 
area of the circle. The solid lines show dispersion relations of 
elementary excitations for the compact support. The broken 
line indicates the chemical potential given by —n 2 t/12. 



Fig. 2. Skew Young diagrams corresponding to the motifs (a) 
01010|1111101|01100 and (b) 01010| 11 11 110| 1 1010. In the case 
(a), there are other possibilities to put numbers or o than the 
one shown in the left figure, while the right figure gives a unique 
combination in the case (b). 



Fig. 3. S(q,w) of the IS supersymmetric t-J model with 16 sites 
and (a) 2 holes and (b) 8 holes. The intensity of each pole 
is proportional to the area of the circle. The solid lines show 
dispersion relations of elementary excitations for the compact 
support. 



Fig. 10. (a) The lslh contribution and (b) the 3slh contribution 
to A(k,u>) of the IS supersymmetric t-J model with 16 sites at 
half-filling. The intensity of each pole is proportional to the 
area of the circle. The solid lines show dispersion relations of 
elementary excitations for the compact support. The broken line 
indicates the chemical potential given by —ir 2 t/12. The meaning 
of black circles is explained in the text. 



Fig. 4. S(q) of the IS supersymmetric i- J model with 16 sites for 
various fillings. 



Fig. 5. S(q, u>) of the NN supersymmetric t-J model with 16 sites 
and (a) 2 holes and (b) 8 holes. The intensity of each pole is 
proportional to the area of the circle. 



Fig. 11. A(k,u)) of the NN supersymmetric t-J model with 16 
sites at half-filling. The intensity of each pole is proportional to 
the area of the circle. 



Fig. 6. N(q,uj) of the IS supersymmetric t- J model with 16 sites 
and (a) 2 holes and (b) 8 holes. The intensity of each pole 
is proportional to the area of the circle. The solid lines show 
dispersion relations of elementary excitations for the compact 
support. 



Fig. 7. N(q) of the IS supersymmetric t-J model with 16 sites for 
various fillings. 



Fig. 12. A(k, ui) of the IS supersymmetric t-J model with 16 sites 
at (a) n = 7/8 (14 electrons) and (b) n = 3/8 (6 electrons). The 
intensity of each pole is proportional to the area of the circle. The 
circles with dark and light shades correspond to A~(k,uj) and 
A~^(k,w), respectively. The solid lines show dispersion relations 
of elementary excitations for the compact support . Th e broken 
line indicates the chemical potential given by eq. (3.23). 



Fig. 8. N(q, uj) of the NN supersymmetric t-J model with 16 sites 
and (a) 2 holes and (b) 8 holes. The intensity of each pole is 
proportional to the area of the circle. 



Fig. 13. A(k,u)) of the NN supersymmetric t-J model with 16 
sites at (a) n = 7/8 (14 electrons) and (b) n = 3/8 (6 electrons). 
The intensity of each pole is proportional to the area of the circle. 
The circles with dark and light shades correspond to A~ (fc, w) 
and A + (k,ui), respectively. 



